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Abstract 

Domain walls, arising from the spontaneous breaking of a discrete symmetry, 
can be coupled to charge carriers. In much the same way as the Witten model 
for superconducting cosmic string, an investigation is made here in the case 
of £7(1) x Z2 — > U(l), where a bosonic charge carrier is directly coupled to 
the wall-forming Higgs field. All internal quantities, such as the energy per 
unit surface and the surface current, are calculated numerically to provide the 
first complete analysis of the internal structure of a surface current-carrying 
domain wall. 
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I. INTRODUCTION 



Domain walls can arise in various grand unified theories (GUT) whenever a discrete 
symmetry is broken by means of a Higgs field. Because they have immediately been shown 
to induce a cosmological catastrophe [0 even if they appear in a very late phase transition, 
their internal structure has not been studied in much details yet, since it was widely believed 
that they could not survive until now. Indeed, with an energy per unit surface of the order 
of the cube of the symmetry breaking energy scale rj say, a single wall crossing the universe, 
even with r\ as low as r\ ~ 100 GeV, would produce an enormous overdensity flw&n ~ 10 8 , 
or, in the case where only small balls were to survive, very large anisotropies in the cosmic 
microwave background radiation would be induced which are not observed 0. Hence, if 
stable walls are to exist in a theory, one must have an inflationary period between the time 
they were formed and now. 

The general belief nowadays concerning domain walls, assuming they were ever produced 
at all, is examplified by the Peccei-Quinn phase transition ||, whose cosmological relevance 
notably for the axion problem is still the subject of open discussion ||]. The idea is that 
even though walls could have been formed, the corresponding phase transition would have 
been preceded by a string forming transition in such a way that domain walls could only 
be bounded by strings. In such a framework, all walls would have had a finite size, huge 
surface tension, and would have evaporated in less than a Hubble time, thereby effectively 
solving the problem. It could therefore appear that studying their internal structure is 
indeed pointless. 

However, just like in the case of cosmic strings, the situation could be rather different if 
domain walls were to have the ability to carry some sort of charge. In the case of strings, 
a current has the immediate effect of breaking the Lorentz symmetry along the worldsheet, 
so that one can consider rotating loops (called vortons || because of their particle-like 
properties, or rings 0). The point is that cosmic strings are believed to scale (see Ref. @ 
for a recent review) because the network of string is dominated by the loops, who eventually 
gravitationally radiate all their energy away. When a current is present, these loops might 
reach equilibrium configurations |5|0 whose classical stability was recently discussed J7],|S[] 
with the result that if no quantum instability exists, the scaling is spoiled and they could 
easily overfill the universe unless they were produced at a very low energy scale (estimated 
at ~ 10 TeV). 

Now if the strings bounding the walls were superconducting, the problem could in fact 
be rather similar, the presence of a domain wall modifying the equilibrium configuration in 
an unknown way, while presumably modifying the constraint. This issue, which can, and 
should be analysed in the framework of Carter's formalism [0] for describing p-branes, is 
still a completely open subject. Another difficulty can arise in the case where strings are 
not current-carrying, but if the wall itself is. Indeed, the point is, as before for the cosmic 
string scenario, that the breaking of the Lorentz symmetry along the worldsheet, whatever 
its intrinsic dimension, allows a definition of rotation, and eventually the recognition of the 
existence of centrifugally supported states. Of course, it is not clear yet whether these objects 
could be formed and reach stable states at all, and therefore their cosmological relevance has 
not been established. However, in the purpose of studying these frisbee-like configurations, 
it is necessary that one knows the relevant internal quantities such as the energy per unit 
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area and the surface tensions: they are explicitely calculated in the present article. 

It may seem that coupling charged (or hypercharged) particles to a domain wall forming 
Higgs fields is a bit arbitrary, but in view of the fact that most topological defects are 
predicted to form in various GUT models where the number of degrees of freedom, including 
scalar, vector and fermion fields is huge, and where the couplings are almost unrestricted, it 
seems fairly plausible. The purpose of this article is thus to present a toy model, similar to the 
Witten bosonic model |KJ for superconducting cosmic strings, where the symmetry breaking 
scheme is simply U(l) x Z 2 — > U(l). This model, much like the Witten model, is expected 
to yield qualitatively relevant results. The work is arranged as follows: after presenting 
the actual model in a first section, we investigate the microscopic structure of such a wall 
and end up by dealing with the abovementionned integrated internal quantities, namely 
the energy per unit area, the surface tensions as well as the surface current. The equation 
of state, relating these quantities, is then computed numerically from the solution of the 
field equations and is shown to share most of the superconducting cosmic string equation of 
state properties |Tl], and in particular the existence of a phase frequency threshold, which 
is discussed in some lenght at the end of the paper. This study of a domain wall model 
thus shed a new light on the general knowledge on current-carrying topological defects by 
showing for instance that a generalisation of the string properties in an arbitrary number 
of dimensions is possible, which in turn give a new understanding of these string features. 
With this idea in mind, we end this article by a derivation of the divergent behaviour of the 
timelike component of the current as a function of the topological defect internal dimension. 



II. WALL MODEL 

Domain walls form whenever a discrete symmetry is spontaneously broken. The sim- 
plest way to achieve that is to break a Z 2 symmetry by means of a scalar whose vacuum 
expectation value shall be taken as (0|y9|0) = ±77, with 77 the energy scale of symmetry 
breaking. This Higgs field may be coupled with hypercharge-carrying fields which we ap- 
proximate by a single complex scalar field S whose vacuum dynamics we require to be 
invariant under some U(l) phase transformation group. In much the same way as was done 
for current-carrying cosmic strings flH||L2]| , we neglect any long range interaction and thus 
assume a global U(l) symmetry thereby emphasizing on the actual dynamics of the 



wall, assuming charge-coupling corrections to be negligible, as was shown to be the case for 



superconducting cosmic strings |T!|. The Lagrangian we shall start with is therefore 



c = - l -\d, v \ 2 - l -\d,n 2 - v^t), (1) 

with the general interaction potential given by 

Vfa S) = - V 2 ) 2 + V - rf) + ^|£| 2 + ^|S| 4 . (2) 

The dynamics given by this Lagrangian include existence of domain walls, i.e., solutions of 
the fields equations that separate domains where (0|y»|0) = +1] from regions where (0|y?|0) = 
—77, and on which therefore (0|^?|0) = 0. For now on, we shall simply write <p for (0|<^|0). 
The wall solution will be a stationnary solution, with the wall locally identified with the 
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(x, y) plane, the various field amplitudes depending only on the third z coordinate. Our 
ansatz is thus 



ip = ip{z) and E = a(z) exp[i(kx — uit)], (3) 

where we have chosen the frame where the spacelike component of the current, defined 
below, is directed along the x direction [this form (|3|) for E can always be attained locally 
by means of a simple rotation in the wall plane]. The conserved current, derived as the 
Noether invariant under phase transformations, is 



J = --^-E* + c.c. = -S d. E*, (4) 
2 58^ 2 ' 



which, with Eq. @ plugged in yields 

Jy, = (kS^ - u5^t)v 2 {z). (5) 
The field equations derived from the Lagrangian (|l|) under the assumptions (||]) read 

§ = t^V - v 2 ) + VAv, (6) 

<]ln [w + 2f( V 2 - V 2 ) + ml + X a a 2 ]a, (7) 



dz 2 

in which we have defined the state parameter 

w = k 2 -cu 2 , (8) 

whose sign reflects the spacelike or timelike character of the current given above by Eq. (|5p 
since 

3^ = wa\z), (9) 

and in the chosen conventions of Eq. (|l|), the Minkowski metric is ^ = Diag { — 1, 1, 1, 1}. 

The possibility of a current in the wall can be seen in two ways. First, one can notice 
that the minimum of the potential, in the actual vacuum, is given by 

(p = ±r] and E = 0, (10) 

and that this minimum is shifted in the wall where ip = to 

\ a m 2 = 2fr ] 2 -ml, (11) 

so a condensate may exist provided 

ml < 2f V 2 . (12) 



Another way to realize that a condensate will in fact appear JT0] in the wall consists in first 
assuming no condensate (E = 0), and solve the perturbative equation for E in the domain 
wall background. For E = 0, the solution of Eq. @ is known: 
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ip = T]ta,nh(-^r]z), (13) 

and setting a perturbation in the form E = a(z)e Mt into Eq. (|7|) yields the one-dimensional 
Shrodinger equation for a 



+ V(z)a = wV, (14) 



where the potential 



V(z) = -2f v 2 [l - tanh 2 (i^^)] + m 2 , (15) 



is negative definite when the condition (|1~2|) holds. Hence, under this condition, £ evolves 
in an attractive potential well, with negative eigenvalues for u 2 . Therefore, there exists 
unstable modes and a condensate forms. 



III. CURRENT QUENCHING AND PHASE FREQUENCY THRESHOLD 

In order to analyse the internal structure of such a current-carrying domain wall, it turns 
out to be convenient to introduce a set of dimensionless functions and variables (, X, Y, w 
and {ai j2 ,3} as 

<p(z)=r l X(C), (16) 

= -^pm (17) 

V -V 



with 

C = vW- (18) 

The state parameter is similarly rescalled into 

W = -^—r. W, (19) 

ml 



and provided we redefine the arbitrary underlying parameters as [11.12 



we get the very simple set of ordinary differential equations 

X" = X[i(X 2 -l) + 2a 2 F 2 ], (21) 

ai Y" = Y[w + 2a 2 (X 2 - 1) + a 3 (Y 2 + 1)], (22) 
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where a prime denotes a derivative with respect to £. 

Two constraints on these parameters arise from the requirement that the theory be 
physically meaningful and consistent with currents flowing along the wall. The condition ([12]) 
for instance, reads in terms of these parameters 

«3 < 2« 2 , (23) 

while demanding that the energy of the wall configuration (if = and X ^ 0) be greater 
than the actual surrounding vacuum configuration (tp = r\ and S = 0) implies 

(a 3 - 2a 2 f < -a 3 . (24) 

The first of these constraints in fact means that there exists a spacelike saturation current 
which cannot be exceeded. To see that this is indeed the case, let us perform an expansion 



of X and Y close to the wall where £ <C 1, in the form fllT 



X~xiC + &C 3 and Y ~ y - a( 2 , (25) 

which satisfy the boundary conditions on the wall worldsheet, and in particular regularity 
of the X field [which accounts for Y'(0) = 0]. Plugging back into Eqs. and fl22"|) gives 

b = Xl (2a 2 y 2 -^), (26) 

and 

a=^[2a 2 -w-a 3 (y 2 + l)}, (27) 

so that because the condensate is actually at its maximum at z — 0, one has a > 0, which 
means 

w < 2a 2 - a 3 . (28) 

Thanks to the requirement (p3|) , we see that the limit applies only in the spacelike current 
case where w > 0, and it reflects the fact that for w = 2a 2 — a 3 , one has y = 0, and 
therefore no condensate, hence no current. So there exist a value of the state parameter 
above which the current quenches to zero. 



On the other hand, investigating the large £ behaviour of Eqs. (|21| ) and (|22| ) yields the 
following asymptotics 

1-X~exp(-C), (29) 

as expected from the knowledge of the true solution (O) in the decoupled case (X a2=Q ~ 
1 - 2e c ), and 



F-expH/^+^C), (30) 
V on 



for positive u; + a 3 , 
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K~cos(J|^% + <5), (31) 
V "i 

for negative w + as, with the special w = — «3 case leading to 

y ~ j^r 1 - (32) 

V "3 

Thus, exactly as in the case of a current carrying cosmic string, there exists a phase frequency 
threshold given by w = —03, or u = m a , above which the integral of the current (|5|) from the 
sheet to infinity diverges. This is therefore not a mechanism depending on the dimension 
of the topological defect under consideration, and can be interpreted as charge carriers 



evaporation from it [fLl |. This phase frequency threshold is discussed more thoroughly at 



the end of the following section where integrated quantities are explicitely calculated. 



IV. MACROSCOPIC QUANTITIES 

For most of the cosmologically relevant calculations with topological defects, it is conve- 
nient to consider them as infinitely thin, and for that purpose, it is necessary to know the 
stress energy tensor and the current as line integrals starting from the wall's worldsheet to 
infinity. For instance, the integrated current reads 



C = 2 J dzyJ\J^\ = 2^H J dza 2 {z) = 1if^fa~ x \v\ J d(Y 2 {(), (33) 



where we have defined v = Sign (w)J\w\ and rescalled it according to Eq. (|1SD; the 
additional factor of 2 is here to account for both sides of the wall. The parameter u, being 
essentially identifiable as k or — u, is readily interpreted and has thus been used as the 
relevant parameter for the plots presented below. 

Another obviously very useful quantity for a macroscopic description of a surface current- 
carrying domain wall is its stress energy tensor 

T r = _ 2g m g upl£_ + g,*c, (34) 

which, in the case under consideration, needs to be diagonalized. It is worth noting at this 
point that even though the existence of a current in the wall indeed breaks the Lorentz 
invariance along the worldsheet, thereby raising the stress-energy tensor's degeneracy, it 
does so through the introduction of one privileged direction. Hence, just like in the string's 
case, there can be only two different eigenvalues, namely the energy per unit area U, and 
the surface tension T. The resulting stress-energy tensor then reads 



(<o) 



( U \ 

-T 

-T 

V 0/ 



(35) 



for a timelike current (for which the spatial isotropy is left unbroken), whereas the spacelike 
current case similarly yields 
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(>0) 



( u 



V 



-u 



-T 



(36) 



0/ 



We shall now calculate explicitely these eigenvalues in the specific case (P under consid- 
eration, and for that purpose, we perform a Lorentz boost in the x— direction in such a way 
that the phase of the current carrier £ reads kz or — ut. In this frame, in which we shall 
for now on remain except when it comes to the lightlike case, one has the energy per unit 
surface 



U = 2j dzT tt = y/v? 3 J d( 



X' 2 + ai Y' 2 + \w\Y 2 + -(X 2 - l) 2 + 2a 2 Y 2 (X 2 - 1) + a s Y 2 (-Y 2 + 1 



(37) 



the surface tension parallel to the current 



Tii = -2 / dzT n 



/ d( 



X' 2 + ai Y' 2 - \w\Y 2 + -{X 2 - I) 2 + 2a 2 Y 2 (X 2 - 1) + a 3 Y 2 (-Y 2 + 1 



(38) 



the surface tension orthogonal to the current 

T ± = -2 j dzTyy = ^V? 3 j d( 



X' 2 + ai Y' 2 + wY 2 + h,X 2 - l) 2 + 2a 2 Y 2 {X 2 - 1) + a 3 Y 2 (W 2 + 1) 

(39) 



while the last integrated component provides a very useful numerical constraint as we shall 
see shortly because 



T z = -2 J dzT zz = ^V/ 3 J d( 



-X' 2 - a x Y' 2 + wY 2 + -(X 2 - l) 2 + 2a 2 Y 2 (X 2 - 1) + a 3 Y 2 (-Y 2 + 1] 

4 2 



(40) 



should in fact vanish identically. This can be checked almost immediately when no conden 

l 

' 2 



sate is present since in that case, one has Xq = tanh^/2, so that X' Q = — \{X 2 — 1) which 



in turn implies 
while the general case gives, with the ansatz @ 



X' 2 + I(X 2 -1) 2 



(41) 



d x T xx = d y T yy = d t T a = 
and finally, conservation of the stress energy tensor d^T^ = yields 

8 7 T ZZ = 0. 



(42) 
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But the boundary conditions one must use are such that asymptotically, the fields take their 
vacuum values, so 



lim a(z) = lim^d z a(z) = ]ha^d z ip(z) = limjy? 2 — rj 2 ) = 0, 



so lim^oo T zz = which, with Eq. (f4"2|) implies T zz = 0. Hence, Eq. ([!(]) provides a 
constraint on the fields, namely 

X' 2 + a x Y' 2 = wY 2 + h,X 2 - I) 2 + 2a 2 Y 2 {X 2 - 1) + a 3 Y 2 {W 2 + 1), (43) 

which is used for numerical purposes since it gives the value of the derivative of X near the 
origin, i.e. x\ with the notation of Eq. (^), as a function of the E field's VEV yo, with 

4 = \ + V>3(1 + 2y 2 ) - 2a 2 - w\. (44) 

Note first that we recover x\ = 1/4 in the noncurrent carrying case, again in agreement 
with the corresponding known analytic solution, and second that Eq. ([43D is not a trivial 
constraint: as numerical integration reveals, the functional U[X(Q,Y(Q] has two extrema 
depending on the field configuration, one of which corresponds to an unphysical maximum, 
whereas the second is indeed a minimum satisfying Eq. (|43|). The numerical program de- 
velopped for solving Eqs. (^Tj) and (|22|) used therefore the constraint ([43D by fixing the 
parameters at the origin with Eq. ([44]) . Two criteria for ensuring the convergence to the 
actual physical solution were thus considered, namely that the solution should be one indeed 
and therefore should extremise U, and the vanishing of T z . 

A last consideration permits an evaluation of the accuracy of the numerical results 
thereby obtained, and it is the final point on the v line calculated for a spacelike current. 
This point corresponds to w — 2a 2 — «3 which, according to Eq. ( p7|) and the discussion 
following it, has no current at all. In that case, all the integrals of Eqs. (37,38,39) are equal 
to U , with 

U = yfctf J d( [X' 2 + i(X 2 - 1)] = 2^ v 3 J dCX' 2 , (45) 
when one takes the solution Xo = tanh£/2, and this is 

U Q = 2^V7 3 / X'dX = -JT^JjX 2 - l)dX = ^^/V/ 3 . (46) 



The condensate therefore must respect 



D ' <\ (47) 



in order to be stable against charge carrier evaporation, with the equality obtained in the 
limit w — > 2a 2 — This in fact also limits the range of variation of w for a timelike current 
for it seems doubtful that a state having U a > U could survive in practice. 

The case of a lightlike current shares with the noncurrent-carrying wall the property 
that the stress energy tensor's eigenvalues are strictly equal. It can usually be set, after 
diagonalization for J M J M ^ 0, as 
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T" v = Uu"u u - T\\x"x v - T ± yVy", (48) 

with the timelike eigenvector {u^u^ = — 1) and x^, the spacelike eigenvectors {x^x^ = 
y^ = 1, and x^ = 0), but for a lightlike current, it reads 

T n» = Uu » u v _ T ^jM _ T ± y»y u - -(u'V + u v x»)u 2 J dza 2 (z), (49) 

where we have set S = a(z) exp[iu)(t — x)]. Upon diagonalization, this reads 

T*" = T\\ {vtv v _ - v%v\ - y^y u ), (50) 

with v± = \{x^ ± u^) the lightlike eigenvectors of T^, and Ty as given by Eq. ( j39f) with 
w = 0. 

Let us investigate more thoroughly the spacelike and timelike cases. The timelike case is 
characterised, as exemplified on Eq. (|35|), by the isotropy of the purely spatial part of T^, 
i.e. Tji =T± = T. As in the string case, one has the Legendre-like relation 

U-T=-vC, z/<0 (51) 

and the now standard formalism developped by Carter |J applies straightforwardly. The 
case of a spacelike current is slightly more involved and perhaps requires more thought for 
each particular cosmologically interesting configuration studied because the spatial isotropy 
of the surface is no longer present since the current picks a privileged spatial direction in the 
worldsheet. However, Eqs. (ftEp, fl57p and ( |59"D show that yet another simplification arises 
from the fact that U = Tj_, i.e. the purely spatial component of the stress energy tensor in 
the direction parallel to the current flow is the energy per unit surface. Setting T = Ty, a 
relation similar to Eq. (O) is obtained in the form 

U — T = uC, v > (52) 

which can be understood in terms of duality between spacelike and timelike currents 0. 
The relevant rescalled integrals are displayed on the figures. 

Fig. 1 represents the energy per unit area and the surface tensions as functions of the 
rescalled state parameter v for a specific set of parameters {a{\ (chosen to yield a generic 
kind of result as well as giving measurable effects), with Fig. l.a showing the variations 
of U(y) and T(p) for a spacelike current- carrying wall having a positive state parameter 
v > 0, while Fig. l.b represents U(y) and T(y) for a timelike current-carrying wall with 
v > 0. Similarly, Figs. 2. a and 2.b show the amplitude of the current (f|) in the magnetic 
and electric regimes respectively. As might have been anticipated, these figure are very 



much like those obtained for a neutral current- carrying cosmic string ||11| , at least in the 
classically stable part of the equation of state, which is definable through the requirement 
that the soundlike perturbation squared velocity c 2 L = —dT/dU be positive. Thus, the 
approximate analytic equation of state proposed in Ref. |13] should be useful also in this 
domain wall context. In fact, the only noticable difference between the wall and the string 
as far as internal structure goes concerns the unstable region: Figs, l.a and 2. a shows that 
almost as soon as the wall becomes unstable with respect to soundlike perturbations along 
the wordlsheet, the wall's stress energy tensor tends to the ordinary wall one, namely the 
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isotropic stress tensor with unique eigenvalue Uq = ^^Jx^r/ 3 . Therefore, most of the current- 
carrying domain wall properties are essentially similar to the string properties. 

Finally, let us remark the following important mathematical property of the surface 
current-carrying domain wall. As is the case for a superconducting cosmic string, it can 
be seen that there exists a phase frequency threshold [TJJ given by w = —m 2 at which 
point the current ( |33|) diverges. For the cosmic string case, the first order pole behaviour 
Cstring ~ (w + m 2 a )~ l was found [[11] whereas the wall case yields Cwaii ~ (w + m 2 a )~ 1 / 2 . This 
is because in both cases, denoting by d the codimension of the topological defect, i.e., 2 for 
a string and 1 for a wall in a 4 dimensionnal background, the current carrier field is seen to 
satisfy Eq. (0) which, far from the topological defect, gives the relation 

A d cr ~ (w + ml)a, (53) 

where stands for the Laplacian in d dimensions: this is simply d 2 /dz 2 in the wall case 
under consideration here, and d 2 /dr 2 + ^—^d/dr in the general case with r the "radial" 
distance to the defect's core. Setting x — kr, with r = z in our wall case and k 2 = w + m 2 
one can extract a as a function of k since for k ^ 0, Eq. (Q) [i.e., Eq. (|53|)1 transforms into 



a- 



dx 2 X dx 
whose solution cannot depend on k 



d 2 a d — Ida , , . „. 

+ ~r~ = a (x), (54) 



The solution to Eq. (54) is well known: 

a ~ Ax'-^Mx), (55) 

with Kq the Bessel function of zeroth order whose asymptotic behaviour is given by Kq(x) ~ 
exp(— x)/y/X- Thus, one finds the general phase frequency threshold behaviour, up to a 
finite part [corresponding to the fact that one has to integrate up to the point where the 
approximation (|53|) becomes valid] 



C oc J r d ~ l dra 2 (kr) 

x h I xK ^ dx 



<x{w + ml)- d/ \ (56) 

with d = 1 for a current- carrying domain wall, d = 2 for a superconducting cosmic string, 
and d = 3 for a charged monopole in a four dimensional background spacetime. It is in 
fact possible to be slightly more precise concerning this divergence: the function k, defined 

asm 

k = 2,-j— = 2, f dV^(xi), (57) 
dw J 

being proportionnal to C, also diverges, and it may be seen that, under the assumption that 

a 3 



Uq ~ ^ near the threshold [see Ref. 0] and Eq. 



11 



k = K f (w) + ALL(w + mlY d l\ (58) 

which is valid for various values of the co dimension d, with Kf(w) the finite part of K and 
A a pure number, calculable in principle by a matching of the asymptotic solution (|55|) 
to the origin and depending on d. Note that the dimension of this function k is given 
straightforwardly once d is known. 
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FIGURES 



Energy per unit area U (solid line) and surface tensions T >0 = T\\ and T <0 = T± (dashed 

lines) as functions of the rescalled state parameter v and in units of ^Jx^rf. Fig. l.a rep- 
resents the equation of state in the timelike case v < 0, while Fig. l.b is for the spacelike 
range v > 0. 

Integrated value of the surface current in units of r] 2 for the same variation ranges as on 
Fig. 1 as a function of v. 
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